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ABSTRACT 


A brief description of how a basic comet ephemeris can be calculated from the orbital elements is 
given. In order to keep mathematics to a minimum, no proofs, derivations or detailed explanations are 
given. 


RESUME 


On décrit une méthode pour calculer l’éphéméride d’une cométe a partir des éléments orbitaux. Afin 
de retenir la mathématique au minimum, aucune épreuve, dérivation ou explication détaillée n’est 
donnée. 


Introduction. In a recent article (Tatum 1982) I described how to calculate the 
position of a comet from measurements made on a photograph. There was one 
little step, however, that I missed out — namely, first find your comet! 

Usually there is no difficulty, because ephemerides for periodic comets are 
published regularly in the Handbook of the British Astronomical Association and 
for new comets in the Circulars and the Minor Planet Circulars of the 
International Astronomical Union. Recently, however, at least partly because of a 
poor mail service, I experienced a long period when the I.A.U. Circulars, for one 
reason or another, were not reaching me, and this at a time when several comets 
were well placed for observation. This made me wonder whether it was possible 
for an ordinary mortal to calculate an ephemeris for himself. 

The answer, I found, was a qualified yes. 

A comet’s orbit is described by certain elements, which are parameters that 
describe its shape, size and orientation. There are two stages to the calculation of 
an orbit and ephemeris. The first is to calculate the elements of the orbit from the 
available observations. The second is to calculate an ephemeris (that is, the 
day-to-day predictions of the right ascension and declination) from the elements. 

The first stage was described by Newton as a problema omnium longe 
difficillimum. Now Newton was no ordinary mortal and I for one am prepared to 
take Sir Isaac’s word for it. The problem of course has long been solved by such 
giants as Laplace and Gauss, but I prefer to leave such heroic calculations to their 
modern successors. 

The calculation of a very precise ephemeris from the elements is no easy 
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problem either. The experts in this field include in their calculations such 
refinements as light-travel time, perturbations by the planets (which are sometimes 
very large) and non-gravitational effects caused presumably by the asymmetric 
evaporation of the material from the comet’s nucleus. 

If, however, one is satisfied with a simple Keplerian orbit without the inclusion 
of such complications, and would like to calculate an ephemeris that is usually 
adequate to find the comet in one’s telescope, the calculation is almost trivial to 
one who possesses a good programmable calculator and knows how to use it. And 
even if the [AU Circulars are being safely delivered, it is still good fun to do it for 
oneself. 

In this article I describe how to do the calculation. I unashamedly give no 
explanation of how to derive the various formulae, since this would result in much 
more mathematics, which would probably deter the very people for whom this 
article is intended. Almost nothing is original, and the material can be found in 
most standard textbooks on celestial mechanics, though I may possibly have added 
to the vast number of methods that have been proposed for the very rapid 
calculation of the true anomaly. 


The Orbital Elements. An orbit can be an ellipse, a parabola or a hyperbola. Since 
a parabola represents a sort of dividing line between an ellipse and a hyperbola, 
and has an eccentricity of exactly unity, one might suppose that parabolic orbits 
never occur. However, most newly-discovered comets have very eccentric orbits 
and are observed over only a very short arc near perihelion. Very often it becomes 
impossible to distinguish between an ellipse and a parabola from such a short arc, 
and an attempt to calculate the major axis of an ellipse is not very reliable. 
Therefore in practice it is very usual to fit a parabola rather than an ellipse to the 
observations of a new comet seen over a short arc. Hyperbolic orbits are very rare. 
An orbit can be hyperbolic if the comet originated outside the solar system and is 
just passing through; or a previously elliptical orbit can be perturbed by Jupiter into 
a hyperbolic one, in which case the comet may be hurled out of the solar system. In 
practice it turns out that the twenty or so hyperbolic orbits that have been observed 
had eccentricities only very slightly in excess of unity, and probably all of these 
were orbits that had been changed by Jupiter from elliptic to hyperbolic. That is not 
to say there are no interstellar comets, but, if there are, they are so few and far 
between that we have not yet knowingly encountered one. In any case, hyperbolic 
orbits are rare, parabolic are common and there are a hundred or so elliptic orbits. 
We describe the size and shape of an elliptic orbit by its semi-major axis a and 
its eccentricity e. In addition an elliptical orbit has a period P associated with it, 
related to a by P” = a? if P is expressed in sidereal years and a in astronomical 
units. All parabolas, of course, have the same shape (e = 1). There is no period, 
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P PLANE OF ECLIPTIC 


Fic. 1—Definition of 1 and w. The perihelion point in this drawing does not seem to be the closest 
point to the Sun, but this is because the ellipse is being viewed obliquely. 


neither is there a major axis, so the size is described by the perihelion distance q. 
The semi-transverse axis of a hyperbola is given the symbol a, and the eccentricity 
is e. There is of course no period. Sometimes the value of a is given as a negative 
number, but in this article I treat a always as a positive quality. 

In addition to size and shape, we must describe the orientation of the orbit. 
Whether it is an ellipse, a parabola or a hyperbola, we do this by means of three 
angles (“Eulerian angles”) as follows: i = inclination of the plane of the orbit to the 
plane of the ecliptic, 0 = ecliptic longitude of the ascending node, w = argument 
of perihelion. 

Figure 1 should help to make this clear. The ecliptic longitude of the ascending 
node is measured eastwards from the first point of Aries ‘VY to the ascending node 
(). The argument of perihelion is measured from the ascending node in the plane of 
the orbit in the direction of the comet’s motion. Both () and w can have any value 
between 0° to 360°. The inclination can have any value from 0° to 180°. A little 
thought will show that negative values of i are not necessary. One can think of 
comets with inclinations greater than 90° as moving with retrograde motion. Other 
conventions for defining i, (2 and w are possible, but the one described seems to be 
the most usual and, I think, the most sensible. 


Auxiliary Quantities Independent of Time. Certain auxiliary quantities are useful 
in computing-an ephemeris. These are denoted here by |, m, n, p, a', b, c, A, B, 
C, a’, B, y. These are functions of the elements i, (2, w and the obliquity of the 
ecliptic 7, but are independent of time. They are useful in transforming 
coordinates from the plane of the orbit to equatorial coordinates. These should be 
calculated as the first stage in constructing an ephemeris, and are defined as 
follows: 
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1 = sin? 0 + cos* 1 cos? i (1) 
m = sin? i (2) 
n= sin 2 71 cos 1) sini cos i (3) 
p = cos Q) cos i (4) 
a’ = (cos* 0 + sin? 0 cos? i)!” (5) 
b = (1 cos* nH + msin*? yn — n)'” (6) 
c = (I sin? n + mcos? yn + n)'” (7) 


The positive square root is to be taken in equations (5) to (7). After this point, the 
quantities /, m, n are no longer required; they were needed only for the calculation 





of b and c. 
Ae eos) Se (- saci (8) 
a a 
pein fn) = cos! [ee @) 
b b 
eeu = sin n) ee — es) (10) 
c C 


It is necessary to calculate A from both parts of equation (8), in order to 
determine the quadrant of A. Both sin”! and cos“! are double-valued, but only one 
value of A, the correct one, will be common to both. A, B and C can conveniently 
be expressed in degrees. After this point the quantities, i, (2 and p are no longer 
required; locations where they may have been stored in a calculator can now be 
freed for other variables. 


Finally, 
a’ =At+o (11) 
B=Bt+o (12) 
y=Ct+o (13) 


and after this point A, B, C and w (as well as 1, m, n, p, i, Q) are no longer 
required. However, a’, b, c, a’, B, y, (which are not functions of time) will be 
needed later. I have used a prime on a’ and a’ so as to avoid confusion with the 
semi-major axis and right ascension respectively. 

Usually we work with the equator and equinox of 1950.0 throughout. In that 
case, n = 23° .44812. 
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Fic. 2—-Definition of v. 


The True Anomaly. The next step in the calculation is to find the true anomaly v at 
time ¢, which is illustrated in figure 2. It is the angle subtended at the Sun between 
the comet at time ¢ and the perihelion point of the comet’s orbit. The comet passes 
perihelion at time 7, so that the time interval between perihelion passage (true 
anomaly = QO) and true anomaly = vis t — T. The procedure for finding v at time ¢ 
is different according to whether the orbit is an ellipse, a parabola or a hyperbola. 


Ellipse. 
First form a quantity M, known as the mean anomaly, from the equation 
2 1(t — 
M = one (14) 
P 
or from 
2(t — T) 
~ aR i 


In equation (14) ¢ — T and P must be expressed in the same units. In equation (15) 
t — T must be expressed in sidereal years and a in astronomical units. The mean 
anomaly is actually an angle expressed in radians. It has a simple geometric 
interpretation, which need not, however, concern us here. 

Now we calculate an angle E, the eccentric anomaly, from the equation 


_ M+ e(sn—E — EcosE), 


E 
l1—ecosE 


(16) 
which is a form of an equation commonly known as Kepler’s equation. The 
eccentric anomaly also has a simple geometric interpretation, though it is not 
essential to describe it here. Some remarks on techniques for the solution of 
Kepler’s equation are given later. 
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cos EF —e 
geese 17 
Cos U <a (17) 


Parabola. 

First form a quantity M defined by 
_ 3nV2(t-—T) _ 13.3286488 (t — 7) 
rr a a 


where t¢ — T is in sidereal years and q is in astronomical units. Then calculate a 
quantity u from the equation 


M (18) 


2u° + M 
Se 19 
304 ne 
Finally calculate the true anomaly from 
v= 2tan! u. (20) 


The sign of v is the same as the sign of M. 


Hyperbola. 

The true anomaly in a hyperbolic orbit can be elegantly expressed in terms of 
hyperbolic functions (sinh, cosh, etc.), but in this article I am not so concerned 
with mathematical elegance as with ease of understanding and with rapidity of 
calculation. Probably some readers are less familiar with hyperbolic functions than 
with the ordinary trigonometric functions (this is certainly true of myself) and few 
calculators have sinh or cosh buttons. The calculation described below may lack 
elegance, but it is very fast and well adapted to electronic hand calculators. 

First form a quantity M from equation (15). Unlike the elliptical case I do not 
know of any particular geometrical significance of M though it is possible that 
there is one. 

Then calculate a quantity u from 


a 2u[e— u(1—|M|—Inyu)] 


21 
u(eu — 2) +e a 
Finally, calculate the true anomaly from 
u(u — 2e)+ 1 
Sy es ee 22 
iid u(eu — 2) +e — 


The sign of the true anomaly is the same as the sign of M. For a hyperbolic orbit the 
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largest possible value of v, approached as t— ~, is cos"! (— e~!). M and u are then 
both infinite. 


Remarks on the Solutions of the Equations. In this section I make some practical 
remarks concerning the solutions of the equations, especially equations (16), (19) 
and (21). It is suggested that, on first reading, this section be omitted. The reader 
can assume that he has succeeded in calculating the true anomaly v and can move 
to the section on heliocentric distance. 

Two small points can be mentioned first. In equation (16) M and E are angles 
and should be expressed in radians. It is perhaps also worth mentioning that in 
equations (16), (17), (21) and (22) the quantity e is the eccentricity, not the base of 
natural logarithms. 

It will doubtless be noticed that equations (16), (19) and (21) at first appear to be 
clumsy and inelegant; by algebraic manipulation they can be re-written in 
apparently much “simpler” forms. For example, equation (19) could be re-written 
as 


u=M-—3w. (23) 


Or again, a u will cancel on each side of equation (21). 

The reader is urged, however great the temptation, not to do this. Considerable 
trouble has been taken to write these equations in a form that leads to the most rapid 
calculation on a programmable electronic hand-calculator. A recasting of the 
equations could increase the amount of calculation ten or even a hundred times. 

In equation (16), E appears on the right-hand side. Likewise u appears on the 
right-hand side of equations (19) and (21). Therefore these equations have to be 
solved by iteration. Convergence is extremely rapid if the equations are not 
“simplified”. To start the iteration, it helps to make a good “first guess”. Graphs to 
supply a good first guess are given in figures 3, 4, and 5. 


The Heliocentric Distance. The heliocentric distance, r, also called the length of 
the radius vector, is shown in figure 2, and is just the distance of the comet from the 
Sun. It is calculated from 


a(1 — e’) 
Ellipse: = —___-— 24 
ee 1 + ecosv a 

2 
Parabola: r= ee. ae (25) 
1+ cos vU 

a(e* — 1) 

Hyperbola: 1c (26) 
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+50° 


+ 30° 
+ 20° 


+ 10° 


(+ if O°<E or M<180°, — if 18O0°<E or M< 360° ) 





= 
2 ee 
o° 20° 40° 60° go° 100° 120° 140° 160° 180° 
360° 340° 320° 300° 280° + 260° 240° 220° 200° 
—E— e 
Fic. 3—Approximate solution of equation (16). 
+0.5 2.0 +5 
+0.4 +1.6 +4 
+0.3 +1.2 +3 
u 
+0.2 +0.8 £2 
+0.1 +0.4 +1 
% +0.4 +0.8 % +4 +8 % + 40 +80 
M 


Fic. 4—Approximate solution of equation (19). 


From this point onwards there is no difference in the calculation of an ephemeris 
between an elliptic, a parabolic or a hyperbolic orbit. 


Geocentric Equatorial Coordinates. We next calculate the geocentric rectangular 
coordinates (X,, y., Z,) of the comet. These are coordinates centred on the Earth, in 
which the xy-plane is the plane of the equator, with the x-axis directed towards the 
first point of Aries, the y-axis directed 90° east of this, and the z-axis towards the 
north pole. To do this we need to transform the coordinates (r, v), which are in 
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Fic. 5—Approximate solution of equation (21). 


the plane of the orbit, first to the plane of the ecliptic using 7, , and w, and then to 
the plane of the equator, using the obliquity n. We have already prepared for this 
part of the calculation by computing the auxiliary angles a’, B, y. 

We also need to know the rectangular equatorial coordinates (x, y, z) of the sun, 
referred to the ecliptic and equinox of 1950.0. These are published for O"UT for 
every day of the year in Section C of the Astronomical Almanac (except for the 
1981 edition), which was known before 1981 as the American Ephemeris and 
Nautical Almanac. 

If we know (x, y, z) the calculation of (x,., y,, Z,) is simple: 


Xe-=ra sn(vut+a')+x (27) 
yo =rbsin(vu + B)+ y (28) 
Z=rcsin(vu+y)+2z (29) 


Finally we calculate the right ascension a, declination 6 and geocentric distance 
A from 


tana = y,/x. (30) 
tan 8 = z,/(x,7 + y,*)'? (31) 
B= ee ye ee) (32) 


The ambiguity in the quadrant of a is determined by inspection of the signs of x, 
and y,, and the calculation is then complete. 
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Example. The elements of the elliptic orbit of Comet Grigg-Skjellerup 1980 II are, according to 
the 1982 Handbook of the British Astronomical Association, as follows: 


T = 1982 May 15.0023 E.T. 
a = 2.958981 A.U. 


e = 0.665683 

i = 21°1366 
Q, = 212°6315 1950.0 
w = 359°3280 


Let us calculate its position for 1982 June 10.0 E.T. (E.T. = Ephemeris Time, which differs at 
present, 1982, from Universal Time by only about 53 seconds.) We recall that, referred to the 
equator and equinox of 1950.0, 


1 = 23°44812. 


I shall proceed with the calculation using the same equation numbers as previously. 
Auxiliary quantities 


| = 0.907781 296 (1) 
m = 0.130027 093 (2) 
n = —0.206 802 170 (3) 
p = —0.785 498 653 (4) 
a’ = 0.980913 661 (5) 
b = 0.995 708051 (6) 
c = 0.215 345924 (7) 
A = 300°847 (8) 
B = 209°791 (9) 
C = 2749862 (10) 
a’ = 300°175 (11) 
8 = 209°119 (12) 
y = 274°190 (13) 


True Anomaly. 
t — T = 25.9977 days. 
But (t — T) in equation (15) must be expressed in sidereal years, and 
1 sidereal year = 365.25636 days. 
Therefore, 
t — T = 0.071 177 sidereal years, 
and therefore 


M = 0.087 862 537 radians (15) 
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Or 
M = 5°0. 
Figure 3 indicates that a good first guess for the solution of equation (16) is 
E = 15° = 0.26 radians. 


We put this value in the right hand side of equation (16) to obtain a better value of E, and we 
repeat this until the iteration converges. In this case it converges to 10 significant figures after two 
iterations: 


E = 149735 923 48 = 0.257 190 383 radians 


E = 149735 63081 = 0.257 185 275 radians. ne) 
Thus we find the true anomaly: 
v = 327199 080 76. 
Heliocentric Distance. 
r= 1.05402 A.U. (24) 
Geocentric Equatorial Coordinates. 
From page C22 of the 1982 Astronomical Almanac we find 
x = +0.204 8638 
y = +0.912 2643 
z = +0.395 5551. 
From equations (27)—(32) we obtain, therefore, 
X- = —0.274 5602 (27) 
ye = —0.008 4583 (28) 
Z, = +0.212 8343 (29) 
a = 12°07706 (30) 
1950.0 
8 = +37°4671 (31) 
A = 0.348 A.U. (32) 
and the calculation is complete. 
Jeremy B. Tatum, 
Department of Physics, 
University of Victoria, 
Victoria, B.C. 
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